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ABSTRACT 

The electromagnetic and weak neutral current matrix elements that enter in the analysis 
of parity- violating quasielastic electron scattering are calculated using a continuum nuclear 
shell model. New approximations to the on-shell relativistic one-body currents and rela- 
tivistic kinematics for use in such models are developed and discussed in detail. Results 
are presented for three closed-shell nuclei of interest: 16 0, 40 Ca and 208 Pb. The current 
work concludes with a study of the sensitivity of the resulting parity-violating asymmetries 
to properties of the nucleon form factors including the possible strangeness content of the 
nucleon. 
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1 Introduction 



Following the discovery of parity violation (PV) in nuclear beta decay, nuclei have played a 
central role in the investigation not only of the properties of the weak interaction but also of 
hidden properties of nucleon and nuclear structure. Prior to the discovery of weak neutral 
currents, investigations of parity violation in electromagnetic nuclear processes were focused on 
the study of asymmetries and polarizations in photonuclear absorption and decays 0. The 
latter were understood as being caused by the parity impurities in nuclear states induced by 
the charge-changing parity-violating part of the nucleon-nucleon interaction mediated by meson 
exchanges @. This parity-mixing of nuclear states continues to be of interest and, for example, 
is at present one of the motivations of proposed experiments to search for nuclear anapole 
moments ||. 

Soon after measurements of neutrino scattering at CERN |4j] implied the existence of weak 
neutral currents, the vector /axial- vector character of these currents was confirmed by the ob- 
servation of parity violation in electron-deuteron scattering at SLAC ||. These experiments 
initiated a new era for the subject of parity violation in nuclei, and since then several ex- 
periments have been undertaken or proposed to measure PV asymmetries in electron-nucleus 
scattering (see Ref. for discussion of the present situation). Two important experiments 
have been carried out so far on complex nuclei, elastic scattering from 12 C at Bates and 
quasielastic (QE) scattering from 9 Be at Mainz ||, yielding PV asymmetries that are, within 
experimental uncertainties, consistent with the Standard Model. Yet there is ample room open 
for speculation, and one of the main interests at present for PV electron-nucleus experiments 
lies in the possibility of extracting information on elusive parts of the nucleon form factors. In 
particular, as discussed at length in Ref. ||, there is at present considerable interest in studying 
various aspects of the nucleon's strangeness content, that is, its strangeness distributions or form 
factors. 

Exploiting the rich variety of nuclei allows one to explore various features of the electroweak 
nucleon current and form factors complementary to those accesible via lepton-proton scattering. 
For instance, the above-mentioned experiment on 12 C performed at Bates J?| was designed to 
filter out the isoscalar electroweak coupling with hadronic matter. Or, as in the experiment 
on 9 Be performed at Mainz ||, measuring parity violation in QE electron scattering has the 
advantage of involving large cross sections, in addition to that of involving various combinations 
of isoscalar, isovector, electric, and magnetic vector (and axial-vector) form factors, which can 
in principle be selected by appropiate choices of nuclei and of kinematics. The price one has 
to pay for this is that the nuclear many-body problem enters the picture, and ambiguities in 
the description of nuclear structure may render ambiguous the extraction of the various single- 
nucleon form factors. Obviously, an unambiguous extraction of nucleon form factors at least 
requires careful evaluation of the underlying mean-field-based description of the nuclear structure 
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involved. 

In Ref. H a study was carried out of the sensitivity to nucleon form factors (particularly 
to axial-vector and strangeness form factors) of PV response functions and asymmetries in QE 
electron scattering from 12 C. The study was carried out within the context of the relativistic 



Fermi gas model (RFG) and then later extended in Refs. JTO], [IT, [12], |13| to more sophisticated 
models that incorporate specific classes of correlations and meson exchange current (MEC) ef- 
fects. In other work reported in Ref. [FJ] an approach arising from models based on relativistic 
mean field theory has been pursued. As discussed in Ref. [[L3|], for appropriately chosen observ- 
ables it appears to be possible to extract useful information from measurements of the QE PV 
asymmetry in nuclei that bears on issues of the strangeness content of the the nucleon and its 
axial- vector, isovector current. Importantly, in the analyses to date the "confusion" from nuclear 
physics uncertainties appears not to obscure such extractions, despite the fact that reservations 
might be raised concerning the nuclear models used flT5| . 

Indeed, the impulse approximation has been shown to be a good approximation for QE 
kinematics and final-state interactions are not expected to affect the PV asymmetries for high 
momentum transfers, so that the models used previously are quite likely to be capable of yielding 
the critical observables for extracting information on the single-nucleon form factors from the 
integrated PV asymmetry. However, it is still important to explore alternative approaches to 
the nuclear dynamics. Only if the appropriate observables are stable with respect to making 
different choices for the nuclear models is it then reasonable to expect that the single-nucleon 
content can be extracted; if different nuclear models yield significantly different results, then 
one must doubt the entire procedure at present, since each model has its special merit and 
no one approach includes all aspects of the problem (relativistic effects, many-body currents, 
correlations of all types — Hartree-Fock, RPA, et cetera). 

Accordingly, in this work we extend the previous work by exploring the various PV observ- 
ables within the context of a continuum shell model (CSM) description for comparison with 
the other approaches mentioned above. In particular, in this work the CSM is applied to three 
closed-shell nuclei, 16 0, 40 Ca and 208 Pb where the approximate expressions used for the currents 
are known to hold best (see Sect. 2). This model has been thoroughly tested by Amaro et al. 



16| , [17] , [18|1 for the computation of (parity-conserving) longitudinal and transverse responses in 



the region of the QE peak. In Ref. |16[] a comparison was presented with results obtained using 
the continuum random phase approximation (RPA) model of Ref. Jl9|] , which has as a residual 
interaction the effective finite-range interaction derived from the nuclear matter polarization 



potential of Ref. [p0|| . It was found [16 that the CSM approach produces similar results at the 
q- values of concern here (the results are practically undistinguishable at q — 550MeV/c). In 
this past work the effects of meson-exchange currents were also investigated and it was seen in 



Ref. |Tj]] that the net effect of MEC is very small except in the very high energy tail of the QE 
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response where contributions of the two-particle emission channel dominate the R T response 
(as an example, for q = 550MeV/c this starts to be important only at u ~ 350 MeV). Thus 
the CSM is expected to be reliable as well for the computation of PV response functions in the 
vicinity of the QE peak. The issue of relevance for the present work is how different the various 
PV observables are in the CSM compared with those that result when other models such as the 
ones mentioned above are employed. Note that previous CSM modeling has been undertaken 
only at the non-relativistic level and only for the case of the EM responses JT7L |TB|; here we 



include relativity in the calculation and consider also the PV responses and asymmetries. 

This paper is organized as follows. In Sect. 2 we very briefly review the formalism for PV 
electron-nucleus scattering reactions based on the material presented in a recent review article 
on the subject ||. The discussion centers around several issues: (1) We review the essential 
expressions for the general nuclear response functions and PV asymmetry; (2) We consider the 
electroweak hadronic currents for on-shell nucleons and drawing on exact answers and various 
approximations to these quantities (whose details are presented in Appendix A) we develop 
the single-nucleon currents and their one-body consequences; (3) We summarize the continuum 
shell model that forms the basis of the present work (again, with some specifics relegated to 
Appendix B and other connections to the RFG model placed in Appendix C). In Sect. 3 we 
present the results obtained for the response functions and asymmetries for the cases of the 
three closed-shell nuclei 16 0, 40 Ca and 208 Pb. The discussions in that section open with a re- 
examination of the approximations made in the present work for relativistic aspects of currents 
and kinematics. In particular, we apply these to the Fermi gas model, where we know the exact 
answer {i.e., the RFG), as motivation for taking the same approach in the CSM. Following 
this the PV asymmetry is explored with the aim of quantifying the nuclear model uncertainties 
inherent in attempts to extract information about the nucleon's strangeness and/or axial- vector 
content from PV QE electron scattering. Finally in Sect. 4 we summarize our conclusions. 



2 Summary of theory 

2.1 Asymmetry and Nuclear Responses 

The general formalism for parity-violating electron scattering has been presented in detail in 
previous work (see Refs. [0, 21 1 and references therein). Therefore, in this subsection we sum- 
marize only those aspects of the reaction that are needed for the discussions that follow. As in 
Ref. || , we limit our attention to the plane- wave Born approximation (PWBA) with single pho- 
ton or Z° exchange. The processes we consider are represented in Fig. 1. Here, an electron with 
four-momentum = (e, k) and helicity h is scattered through an angle 9 e to four-momentum 
K ffl = (e', k'), exchanging a photon or a Z°, where the transferred four-momentum in the pro- 
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cess is given by = (K — K')^ = (lu, q). We follow the conventions of Ref. [|22|, including the 
invariant phase space factors that are neglected in some other work. 

The basic quantity of interest here is the helicity asymmetry, defined as the ratio involving 
the difference and the sum of the electron scattering cross sections for positive and negative 
incoming electron helicities. The difference is parity-violating, while the sum is twice the usual 
parity-conserving electromagnetic (EM) cross section. When considering single-arm scattering 
of longitudinally polarized electrons from unpolarized nucleons and nuclei, the leading-order 
PV contribution arises from interferences between the two processes shown diagrammatically in 
Fig. 1. PV effects can also come from Z° exchange alone, but these are much smaller than the 
interference terms considered and thus are neglected here. Note that the 7-exchange amplitude 
(a) is purely vector, whereas the Z°-exchange weak neutral current amplitude (b) has both 
vector (V) and axial-vector (A) contributions. The helicity-difference asymmetry may then be 
written 

da + da~ 

A= MMZMdl = Ao ^L } (!) 
da + da W em 

dVL'de' + dVL'de' 

where the scale is set by the term 

A - ° m (2) 

with G the Fermi coupling and a the fine structure constant. 

The nuclear physics (hadronic) content in the problem is contained in the ratio of the parity- 
violating responses (W pv ) to the familiar parity-conserving EM responses (W em ). Within the 
PWBA, these are given by 

Wem = V L R L + V T R T (3) 

Wpv = v l R l av + v t R t av + vt>RVai (4) 

where L and T denote longitudinal and transverse projections with respect to q, respectively. 
The whole dependence on the electron kinematics is contained in the usual kinematic factors, 



v L = Qy q \ VT = tan 2 (fl e /2)-Q7(2g 2 ) andv T , = tan(0 c /2) ^Jtan 2 (9 e /2) - Q 2 /q 2 . The subscript 
AV in the PV responses denotes interferences of axial-vector leptonic currents with vector 
hadronic currents and the reverse when written V A. 

The five hadronic responses in Eqs. can be obtained from the components of the EM 
hadronic tensor W££ in the case of the two pure EM responses R L and R T , and from the hadronic 
tensor W^i nc which involves interferences of EM and weak neutral currents (NC) for the three 
PV responses, R\ v , R^v an d Rva- By considering the laboratory coordinate system with the 
z-axis along q, the y-axis along kxk' and the x-axis in the electron scattering plane, the various 
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hadronic response functions may be written as follows, 

R L = wZ = EW/-^ - u,)\(f\p™\i)f (5) 
r t = w:z + wz = EE*(%-^-«)K/ra»*>i a . (6) 

where the sum runs over all the final unobserved states and an average over initial states is 
performed. The terms p em and J e ^ represent the Fourier transforms of the EM nuclear charge 
and transverse current operators, respectively. Current conservation has been used to eliminate 
the z-component of the hadronic current. 
For the PV responses we can write 

Rav = -y^e °/„c = -fEE^/-^-^)(/IP em K)*(/l(P nc VK) (7) 



dT _ y A ( \xrxx | Txryy 
U AV — 2 y V em/nc + VV cm/nc 

= -fEW ~ E i- ")(f\JTW ■ (f\ (J Q ±)v \i) 

f i 



R t v'a = - l -f(wz /nc -wz 



/nc 



, E ES(Ef -Ek-uj) [(f\J cm \i)* x (f\ ( J nc ) A |<>] • * . (9) 
/ / > q 

Here the indices V and A refer to the vector and axial-vector contributions to the hadronic 
neutral current. The terms {p nc ) v and (J r j°)y are the Fourier transforms of the NC vector nuclear 
charge and transverse current operators, respectively, while (J nc ) A is the Fourier transform of 
the axial-vector contribution in the NC operator. 

The nature of the leptonic vertex appears here in the vector and axial-vector electron cou- 
plings, gv and qa- The Standard model for the electroweak interaction at tree level is assumed 
in this work, 

gv = -1 + 4sin 2 ^ ^ -0.092 (10) 
9a = 1, (11) 

and 9w the Weinberg angle (sin 2 6w — 0.227, as discussed in Ref. 0). 
2.2 Electroweak Hadronic Currents 

As noted from Eqs. (fjH(|5|), in order to calculate the EM and PV responses one needs to 
evaluate the pure EM and weak NC hadronic currents. Following the detailed study presented 



in Refs. ||, |2lJ and eliminating terms involving (c, 6, t) quarks, whose contributions to nuclear 
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matrix elements of {J™)v/a are supressed, the nuclear vector and axial-vector NC operators 
may be written as follows, 



Z ( y- l) j° m (T = i) + V3&- 0) j; m (T 

ei T=1) 4 3) +d T=0) 4 8) +d 0) 4 s) > 



(12) 
(13) 



where Jf™ 



(T = l)/J° m (T = 0) are the isovector/isoscalar EM currents, and = sj^s is 
the strange quark contribution in the vector NC current. The operators Ajf* are given by 
AW = gA a 7 A( 75g/2, where g represents the triplet of quarks (it, d, s), and the A a , a = 1...8 
are the Gell-Mann SU(3) matrices. The strangeness content in the axial-vector NC current is 
given by A^} = s^^s. Finally, the £y's and ^'s are couplings determined by the underlying 
electroweak gauge theory (see Refs. || |2l|). We use the minimal Standard Model tree level 
couplings 

1 £i' 1 



AO) 



(T=0) 



(T=l) 
V 



-4 sin 2 9 W 
- 4 sin 2 9 W 



AO) 

<,A 
JT=0) 







(14) 



AT=1) 



Single-nucleon matrix elements of the EM and weak NC currents shown previously are re- 
stricted by Lorentz covariance, together with parity and time reversal invariance to the following 
forms, 

F 2 



(N(P')\J™\N(P)) = u(P') 
(N(P')\{J* C ) V \N(P)} = u(P') 
(N(P')\(j?) A \N(p)) = u(P') 



Fi7m + i 



F 2 



GaIv 



2M 
M 



u(P) 
u(P) 



(15) 
(16) 

(17) 



where u(P) and u(P') are the single-nucleon wave functions properly normalized, Q = P' — P 
is the four momentum transfer to the nucleon and M is the nucleon mass. From Eqs. (|l~2"l,|l3"l) 
and (|T5D-(fPTD one can write, 

Fa = ^F^rs + VS^F^ + ^F^, a = 1,2 (18) 
G a = ^GWrs + er^GW+^GW, a = A,P, (19) 

where F^ T=0 '^ denote the isoscalar and isovector EM Dirac and Pauli form factors of the nucleon, 
the G^ 3,8 -* are the triplet and octet axial-vector form factors, and F^ and G^> are the vector 
and axial-vector strange-quark form factors. In Eqs. (f0|) and fli~9|) the terms involving t% are 
isovector while the rest are isoscalar. We will mainly use throughout this paper the Sachs form 
factors defined as: Ge = Fi — tF 2 and Gm — Fi + F 2 , with r = \Q 2 \/4M 2 . Analogously one 
can also define Ge and Gm from F\ and F 2 . Also, as discussed in Appendix A and Ref. |J, the 
pseudoscalar contributions are absent in PV electron scattering. 
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In this work we follow the usual procedures of employing the on-shell single-nucleon currents 
to construct one-body EM and NC current operators. However, in contrast to some past work, 
we devise operators that retain important aspects of relativity. In Appendix A expressions 
for the exact on-shell operators for use between two-component spin spinors are summarized 
together with the approximations made in the present work when using these operators within 
the CSM. Below we repeat these approximate expressions for use in building the one-body 
nuclear current operators needed in treating PV electron scattering; specifically, only the vector 
time projections (EM and NC) and transverse vector (EM and NC) and axial-vector (NC) 
projections are required. As also discussed in Appendix A, since we consider only doubly- 
closed-shell nuclei, it is also a reasonable approximation to drop the spin-orbit contributions to 
the time projections, as we do in the following. With dimensionless variables k = q/2M and r 
defined above we then have 



(P nc )v 

jrem 

(J nc )v 
(J nc ) A 




r p + p ' ° x q 

<j£ - „ , r %(jm- 



Ge 



2M 
P' 



2M 
VI + tG A (t, 



+ iG 




(20) 
(21) 

(22) 

(23) 
(24) 



where for the three-vector components of the current only the transverse projections are to 
be employed — the longitudinal projections of the vector currents are given by the continuity 
equation and no longitudinal projections of the axial- vector current are required in descriptions 
of PV electron scattering. Therefore, within the IA, the hadronic (one-body) current operators 
are given by the following: 



Time Components 



p cm = 



(P nc )v = 



K A 
K A 

_ V e iq ' r fc 



1 +T. 



Gl- 



+G'e- 



T k ~ 1 

^+G r i 



To 



'E~ 



(25) 
(26) 



Spatial Components (transverse projections only) 



jcm = 



k k 2M k 

l + r* 



+ G 



M~ 



+ \G E ^ + G E - 



[2V k + iq) 



q x a 



(27) 
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IV 



„t A e iq ' rfe 
-i^Y- 



1 + r; 



G 



1"^ 



M 



q x cr 



+ G p P 



^ 1 + r ; fc 
2 

E 

fe=i 



3 +G" F ^i 



(2V fc + iq) 



iqr fe 



r .4" 



~ 1 — T fc 

<-J 4 ~ r <~r^ 



:28i 



(29) 



where is the position of the fc-th nucleon, while cr k and Tg are the Pauli spin and (third 
component) isospin matrices for the fc-th nucleon. As is standard we refer to the contributions 
proportional to Gm and Gm as magnetization contributions, and to those proportional to Ge 
and Ge as convection contributions. 



2.3 Nucleon Form Factors 

As seen in the previous section, the EM and NC current operators are given in terms of the 
familiar EM form factors, G^™ and G^ n , and the weak NC ones, G^ n , G V m and G V A . Note, 



however, that by using Eqs. (18,19) one can express these weak NC form factors in terms of the 
purely EM ones G^ n M , the axial- vector cases, G p A ,n , and the form factors G^ MA that enter when 
the nucleon has nonzero ss strangeness components. At tree level we obtain for the nucleon 
weak form factors the following expressions (see Ref. ||): 

G P e™m = (1-4 sin 2 8w)G p £'m ~ G e ,P m ~ Ge,m (30) 
&f = TlGf+G { i. (31) 

From these results one easily sees that the strangeness form factors have an isoscalar charac- 
ter. We must also have G^ = for Q 2 = u 2 — q 2 = 0, since the nucleon has no net strangeness, 
whereas the r = values of the remaining strangeness form factors may be nonzero at Q 2 = 0. 
In the case of the axial-vector weak neutral current form factor G A , all of the isoscalar depen- 
dence is contained in the strangeness form factor G A , while the remaining term G A is purely 
isovector (T = 1). Note that this result is only valid at lowest order in the Standard Model 
where we have ^~°' ) = 0. Therefore, in the absence of strange quarks, the isoscalar axial-vector 
form factor G A ~ would vanish at tree-level (for a more detailed discussion of this subject see 
Refs. |], [?]]]). There exist different conventions for the NC couplings; the inter-relationships 
between them are presented in Ref. |J. 

In order to evaluate potential PV electron scattering experiments, it is mandatory to charac- 
terize the Q 2 -dependence of the various nucleon form factors. In Sect. 3.3 we present a detailed 
study of the sensitivity shown by the asymmetries to specific aspects of the nucleon form factors. 
In particular, we compare the results obtained by calculating the EM form factors using the 
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so-called Galster parameterization (see Refs. 



my. 



LTR — <_TT/ 



p v G 



D 



G 



M 



PnTGy^ n 



(32) 
(33) 
(34) 
(35) 



with G$(t) = [1 + A^r]" 2 , £ n (r) = [1 + A n r] _1 , \y = 4.97 and A n = 5.6. In Sect. 3.3 we consider 
variations of A n . 



Following Refs. |g, |21| the isovector axial-vector form factor is given by 

—QaGaj 



(36) 



with qa = 1.262, and G A (r) = [1 + A^t] -2 is the axial-vector dipole form factor. In this work 
we have assumed the value \® = 3.32 corresponding to a dipole mass of 1 GeV. Effects due to a 
renormalization of can be quite significant in neutral current processes. In particular, in the 
case of electron scattering, uncertainties in theoretical estimates of higher-order contributions 
may be on the order of ±10% (for a more complete discussion of these issues see Refs. || |2l|). 
Again, in Sect. 3.3 we return to consider variations in the magnitude of gA- 

Finally, generalizing the Galster parameterization, the (purely isoscalar) strange form factors 
can be written as 



G w 



G 



M 

w 



PsTGyt^E 



(37) 
(38) 
(39) 



The values chosen for p s and p s will be discussed in Sect. 3.3, whereas the remaining form factor 



Gjl will be ignored in the present work as its effect in QE electron scattering is expected to be 



very small ||. One should note that the choice given by Eqs. (|37|)-(|39|) is somewhat arbitrary, 
except for the r — > behaviour of G^\t) dictated by the fact that the nucleon has no net 
strangeness. Given that no experimental information at all exists on G^ M , while very little 
is known about G A (see Ref. ||), Eqs. (|37f )-([3~9"l) simply show a straightforward extension of 
the pure EM form factor equations to characterize the present lack of knowledge in these form 
factors. Results of future measurements may require a different choice of parameterization. In 
this work, we describe the leading, non-trivial Q 2 -dependence of G^ea the parameters p s , 
p s and g A s \ respectively, and possible deviations of the high-r dipole fall-off by the functions 
^!ema t na t show a r-dependence in the form £$ma{ t ) = (1 + ^ema t )~ 1 - The parameters 
X~em A are ' as y e ^' unconstrained. Our focus in the present work is entirely on the region where 
r < 0.3. 
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2.4 Nuclear Shell model 



As has been already mentioned in previous sections, in this paper we restrict ourselves to the 
Impulse Approximation (IA), i.e., we only consider one-body current operators. The study 
of MEC will be presented in a forthcoming publication. In this context, apart from different 
prescriptions for the nucleon form factors, one must also model the nuclear structure in obtaining 
the nuclear matrix elements that enter in the various EM and PV response functions. In this 
subsection we briefly discuss the nuclear structure model used to describe the ground and excited 
states of the nuclei under consideration. In this work we restrict our attention to nuclei with 
doubly-closed shells, 16 0, 40 Ca and 208 Pb, taking the initial states \i) = |0) to be given simply 
by the degenerate ground-state wave functions of such nuclei, whereas the final states |/) are 
represented by lp-lh excitations (within the IA). 

In the following we summarize the basic features of the CSM. The particle and hole, pro- 
ton and neutron wave functions are solutions of the radial Schrodinger equation with a mean 
potential of Woods-Saxon type 

™ - + ( 1+exP( ;- fl „)/ a J « • ' + c*» 

where m n is the pion mass, and Vcoui is the Coulomb potential of a uniform charge distribution 
with radius -Rcoui and charge (Z — l)e. Obviously, the Coulomb potential only interacts in the 
case of protons. The parameters are fixed to have a good description of the empirical single- 
particle energies near the Fermi level and of the experimental ground-state charge density. We 
take -Rcoui = Rq- The rest of parameters are given in Table 1 and have been taken from 



Refs. [[18|, 7A 1 In previous work ITB, O, TE\ we have analyzed the sensitivity of the Rl and 



Rt responses to variations of these parameters. 

The asymptotic form of the particle radial wave functions of momentum k = \JlMe and 
angular momenta l,j is 



Rdj(r) ~ \ ^—-sm(kr-r]log2kr-l^- + ai + 8i j ), (41) 
V ixk r Z 

where Sij is the nuclear phase shift, o\ is the Coulomb phase shift (which is absent for neu- 
trons), and r) is the Coulomb parameter MiZ — l)e 2 /k for protons and zero for neutrons; the 
normalization of the wave functions is then 

(R dj \R enj ) = 6(e-e'). (42) 

Special care is needed when the energy is small or very large, when the angular momentum I 
is very large, and when the charge Z is very large, since then the computation of the Coulomb 
wave functions is known to be somewhat delicate. 
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In this approach, since we only consider inclusive electron scattering, the nuclear final state 
is taken to be an eigenstate of the total angular momentum of the many-particle system, i.e., 



\f) = \a,J) = \p,h~ 1 ;J) 



|0>, 



(43) 



where aj, creates a particle in the continuum, namely state \p) = \e p l p j p ) with e p > 0, while 
b\ creates a hole in the core |0) corresponding to destroying particle state \rihlhjh)- Then, for 
a fixed value of the excitation energy the nuclear excited state is calculated by summing over 
all of the particle- hole pairs (J2 a ) that are allowed by the angular momentum selection rules 
and energy conservation. The calculation of the responses requires the sum over all the excited 



states (see Refs. |L6], [L7| [18], ^6) for details concerning the calculations). 

When considering initial and final nuclear states with good angular momentum quantum 
numbers it is convenient to expand the hadronic currents in terms of multipole projections of 
the charge and transverse three-vector current operators: 



rem 

' ' tin 



J>0 

J>1 



If J + mX, 



MJ 



(44) 
(45) 



with m = ±1. We use the notation J = \/2J + 1. Note that the longitudinal component of the 
three-current can be obtained from the charge current by using current conservation. 

The above multipole decompositions can be applied to the EM currents, p em , J cm , as well 



as to the vector and axial-vector weak neutral currents, (p 11 



IV; 



The 



corresponding Coulomb and transverse electric and magnetic multipole operators are given as 
usual by 



Ml 



rpEJ 

m 



iMJ 



drjj(qr)Y r j n (r)p(r) 

- [dr{Vx[jj( q r)Y^(r)]}.J(r) 
a J 

drjj{qr)Y™j{r)-J{r), 



(46) 
(47) 
(48) 



with jj(qr) the spherical Bessel functions, YJ"(f) the spherical harmonics and l^j(f ) the vector 
spherical harmonics. Here p(r) and J(r) stand for any of the nuclear EM, vector NC and axial- 
vector NC current operators. 

Introducing the above multipole analysis in the general expressions for the response functions 
(Eqs. (|)-(H)) an d making use of the Wigner-Eckart theorem, it is straightforward to obtain the 
following expressions: 



R L (q,uj) = 47r££5(£a-£o-a;)|(a, J\\M t 

J>0 a 



(49) 
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R T {q,u) = 47rJ2J:8(E a -E -u)\(a,J\\T^ + iTM J p}\ 2 (50) 

J>1 a 

R L AV (q,u) = -2ng A ^^5(E a -E -u)(a,J\\MJjO)*(a,J\\(Mi) ||0) (51) 

J>0 a 

R T AV (q,uj) = -2ng A J2J2S(E a -E Q -u)[(a,J\\T^\\0r(a,J\\(T^) \\0) 

J>1 a 

+ (^J\\T e A l J \W(a,J\\(T^ J ) v \\0)} (52) 
JJ&faw) = -2ng v J2^8(E a -E Q -u J )[(a,J\\T^\\or(a,J\\(T^ J ) A ^^^^ 



,/>l a 



+ 



«^l|T c A m J ||0)*(«,J||(T n ^)j|0)], (53) 



where we have used the notation, M^ m , and T^ J for the EM multipole operators, , 
(T^ c jS j v and (T^ J ^J v for the multipole operators associated with the vector NC four-current, 
and (T^ C J ^J and (T^ J ^J for the ones corresponding to the transverse axial- vector NC current. 

Note that the Coulomb and electric operators have parity (— 1) J in the case of EM and weak 
NC vector currents, but they have parity (— 1) J+1 for the axial- vector current. On the other 
hand, the magnetic operator has parity (— 1) J+1 for the EM and vector NC currents, and (— 1) J 
for the axial- vector one. 

Explicit expressions for the reduced matrix elements appearing in Eqs. (|49| |53|) are given in 
Appendix B. It is worth pointing out that the following relationships hold for the approximate 
current matrix elements: 



(a 



, J\\ (T* J ) JO) = Jl + l^{a,J\\TZ J \\0U, 



t G M 

T Gm 



«^ll(^ J )j|0) m a g (54) 



(a,J\\(T™ J ) A \\0) = Jl + ^KJIIT^IIO), 



T G M 

tG 



(a,J||(T n ^)j|0) mag , (55) 



M 

where here "mag" denotes the magnetization current. Theoretically the sums in Eqs. (J4Q|) — (p3|) 
extend to J max = oo, while in practice we add together only a finite number of multipoles. In 
general we use a value of J ma x that is greater than that needed for convergence, J CO nv, specifically 
the values given in Table 2. 

We must take into account the resonant part of the excitation spectrum at low energy. 
In addition we need to integrate the responses to compute the sum-rule asymmetry given in 
Eq. (|57D . Thus we need to compute the responses with a small step size, from emission threshold 
up to a sufficiently high energy uj max that the responses become negligible compared with their 
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peak values. In particular for q = 300, 500 and 700 MeV/c we use as step sizes 1, 2 and 5 MeV 
and u max = 200, 350 and 500 MeV, respectively. 

At low u resonances in the Woods-Saxon potential produce rapid oscillations and sharp peaks 
in the responses, whose overall contribution we take into account by convoluting the responses 
with a Gaussian weight function f(w) oc exp(— u 2 /T 2 ) of width T (in all cases we have used 
T = 5 MeV). So what we show in the plots is not actually the shell model response Rsm(q, w), 
but the "smoothed" response Rr(q,ou), given by 

/CO 
duj'f{uj -L}')R SM (q,u'). (56) 
-oo 

In the next section we compare the results obtained within the present model to those 
obtained with the RFC To facilitate the discussion there, we summarize in Appendix C the 
expressions for the five response functions within the relativistic Fermi Gas Model together with 
the non-relativistic limits. 

3 Results 

In this section we present and study in detail the results obtained for the response functions and 
asymmetry defined in previous section and appendices of the present paper. An important goal 
in this work has been to explore the sensitivities that the various observables display to different 
nuclear models, as well as to properties of the single-nucleon form factors. Concerning the former, 
our motivation has been to investigate whether or not shell effects show up in the observables 
under consideration and to explore these nuclear model dependences using an approach that 
differs from that employed in past work in this area. In particular, as discussed above, we have 
argued that the CSM augmented by the incorporation of an approximate relativistic one-body 
current and relativistic kinematics should provide a reasonable model for inclusive medium- 
energy QE electron scattering. Our aim in the present section is to present results that support 
this approach and then to proceed to investigate aspects of the single-nucleon structure physics 
as revealed in PV QE electron scattering. 

To set a clear scale for the model dependences (and therefore the ambiguities) in the re- 
sponse functions and PV asymmetry, we have performed calculations for different values of the 
momentum transfer ranging from q = 300 MeV/c up to 1000 MeV/c using the CSM described 
in Sect. 2.4. In so-doing we deal with situations in which non-relativistic calculations are ex- 
pected to be reliable, and others in which a relativistic analysis of the reaction mechanism is 
mandatory. Results for three closed-shell nuclei, 16 0, 40 Ca and 208 Pb, are presented. Naturally, 
the PWBA used in this work is not directly applicable in modeling the response functions for 
nuclei as heavy as lead. In such high-Z cases the electron distortion can be taken into account 
by obtaining numerical solutions to the Dirac equation with a Coulomb potential. In the present 
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paper our main focus is placed on the PV asymmetry which, being proportional to a ratio of 
helicity difference and sum cross sections, is less affected by distortion. Accordingly, for the 
model-to-model comparisons made here the PWBA should prove adequate. 

Throughout this work, the nucleon's EM form factors are taken to be given by the simple 
parameterizations discussed in Sect. 2.3. Since it is expected that these functions will be studied 
intensively at most medium- to high-energy electron scattering facilities over the next few years, 
we have not pursued any further analysis of the sensitivity of the PV asymmetry to this aspect 
of the problem other than to present one figure showing the sensitivity to variations of G\] this 
issue can easily be re-examined when new data become available from measurements of parity- 
conserving polarization-dependent electron scattering from the proton, deuteron and 3 He. 

We begin by reproducing the observables in the RFG model studied in previous work ||, 
lPj p| , |12| , |i"3l (see Appendix C) to quantify further the approximations made in our present 



treatment of the CSM. Importantly, we show below that the current developed in Appendix A 
provides an excellent approximation to the exact on-shell results when used within the context of 
the Fermi gas model, confirming those past studies. The same approach is then used in the CSM 
to proceed to an investigation of the nuclear model dependences expected for PV observables. 

3.1 Quasielastic responses in the RFG model 

The RFG model has been summarized in Appendix C and several approximations to the single- 
nucleon current discussed in Appendix A. In the present subsection we present results using this 
model to help in assessing the quality of those approximations for use later in the CSM. 

In Fig. 2 we show the EM response functions R L and R T for the case of 40 Ca at q — 500, 700 
and 1000 MeV/c. The solid curves result from using the fully relativistic model obtained with 
the exact on-shell current operators given in Eqs. (|65|-[r0|), i.e., using the formalism presented 
in Ref. ]l(J. Those shown as long-dashed curves correspond to making the approximations 



for the current operators given in Eqs. (|20| ) and (^). Moreover in this latter case we have 
begun with the non-relativistic Fermi gas model for the rest of the problem (i.e., other than the 
currents) and proceeded to make the A = u/2M — > A(l + A) replacement discussed at the end 
of Appendix C, accounting approximately for relativistic kinematics in the energy-conserving 
5-function. As is clearly seen by comparing the solid and long-dashed curves, and in accord with 



past work |TT| , [12], [L3|, when both approximations are made the agreement is excellent, verifying 
that the use of the particular approximate currents and kinematics assumed here yields results 
that differ by only a few percent from the exact RFG results. 

Also shown in Fig. 2 (as short-dashed curves) are responses having the relativized kinematics, 
but now using non-relativistic currents. In particular for these the factors n/y/r and \fr j k in 
Eqs. fl20D and (|22|), respectively, are set to unity. Without the relativistic effects in the currents 
the results for R L (R T ) are too small (large), as expected. We note that conventional treatments 
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of the single- nucleoli current (see, for example, Eqs. (D4) and (D5) of Ref. |57[| ) often take to <C g, 
in which case k = sfr. Moreover, as discussed in Appendix A, the non-relativistic expansions 
are made not only with respect to rj, but also k, A and r. Figure 2 demonstrates the magnitude 
of the error to be expected for QE scattering if one makes such approximations. 

Finally in Fig. 2 are also shown results where the replacement A — > A(l + A) is not made 
(dot-dash curves with the relativistic corrections to the current, dotted curves without). With 
non-relativistic kinematics the QE peak occurs near u = q 2 /2M (A = k 2 = r(l + r)), whereas 
with relativistic kinematics it occurs near uo = \Q 2 \/2M (A = r), namely, at a smaller value. 
As discussed in Ref. ||10|| , the relativistic model yields a smaller width to the QE response than 
does the non-relativistic model. Also, the four-momentum transfer is decreasing as u increases 
for fixed q and consequently the single-nucleon form factors (those that are simply proportional 
to dipole form factors) increase with increasing uo. This is the underlying reason for the large 
enhancements seen in the non-relativistic results shown in the figure. 

The importance of maintaining relativistic kinematics in the energy- conserving ^-function is 
obvious from these results. Indeed it would be a very poor approximation to ignore either of 
the "relativizing" steps made in the present work. Clearly as the momentum transfer reaches 
the medium-energy domain of interest in this work the effects of relativity in currents and 
kinematics (at least) are rather important. Typically one would incur significant error in using 
the completely non-relativistic FG model and thus, by extension, one must doubt any other 
model which does not incorporate the effects discussed in the present work. The CSM model 
discussed below has been developed with this in mind: we use only the relativized current 
operators and always make the replacement u — > u)(l + u/2M) in the energy-conserving 5- 
functions. We now turn to the results obtained using that model. 



3.2 Quasielastic responses in the CSM 

An important aspect in the study of nuclear reactions is to determine as precisely as we can the 
nuclear physics dependences inherent in the observables. On the one hand, as in much of the 
present work, we may be interested in minimizing such dependences with the hope of extracting 
information about the nucleon's various form factors; on the other hand, it may be that certain 
observables contain unusually large sensitivity to some specific issue of interest in studying the 
nuclear many-body problem. The focus in the present work is largely to address the former of 
these two points of view specifically for the PV electron scattering observables. In this regard, 
the study of ratios of responses such as the PV asymmetry is favored, since they usually turn out 
to be less sensitive to the nature of the underlying dynamical assumptions made when compared 
with the response functions (or cross sections) themselves. 

In particular, here we compare the results of calculations using the CSM with well parameters 
given in Table 1 with those of the RFG model for the response functions, the PV asymmetry 
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and the sum-rule asymmetry, the last being defined below following previous studies in this 
area [13|. For the RFG results shown here we use as Fermi momenta the following: pp = 
215 or 225 MeV/c for 16 0, 235 or 245 MeV/c for 40 Ca and 230 or 250 MeV/c for protons in 
208 Pb. In each case these have been chosen as representative for the three nuclei, with a smaller 
value that is favored somewhat when ground-state properties are emphasized and a larger value 
that yields rough agreement when the RFG model predictions are compared with experiment 
(this presumably arises since QE scattering at intermediate energies involves the ground state 
and a rather excited particle- hole state, both of which the model must recognize). For 208 Pb 
we show results using two different assumptions for the relationship between the proton and 
neutron Fermi momenta, (1) with pp n = (N/Z^^ppp = 265 or 288 MeV/c for the two cases, 
corresponding to keeping the volume of the neutron and proton gases equal, and (2) with 
PFn = Pf p = 250 MeV/c, corresponding to equal densities for the two gases. As usual we take 
PFn = PFp for 16 and 40 Ca. In the next subsection we show results for the PV asymmetry that 
indicate the very weak dependence of this observable on the Fermi momenta chosen, confirming 
past work ||. 

The resulting EM and PV response functions are displayed in Figs. 3-8 for q = 300, 500 and 
700 MeV/c (the PV results given here are shown in the absence of strangeness; see the later 
discussion in the next subsection). Specifically, in Fig. 3 we show the EM response functions for 
16 0, followed in Figs. 4-8 with results for the complete set of five electroweak inclusive responses 
for the cases of 40 Ca and 208 Pb. The inter-comparisons of R L and R T for oxygen and calcium 
are typical in that the full sets of responses in these two cases are very similar. Thus, for brevity 
we have focused on calcium and lead to show the nuclear model dependences obtained in going 
from light to heavy nuclei while omitting the PV responses for oxygen. 

From the results in Figs. 3-8 we observe the following: 

1. Overall, the RFG yields a reasonable description of the responses when the momentum 
transfer is high enough. The region where the RFG is Pauli-blocked (the q = 300 MeV/c 
case in the figures presented here) is where the agreement is the worst, as expected, since 
the RFG is known to become a poor approximation in this regime. 

2. On closer examination, one sees that the CSM produces tails extending below and above 
the response region where the RFG is nonzero. Again, this is expected, since the latter 
model has a sharp cutoff to its momentum distribution while the former has all momentum 
components present, these being governed by the initial- and final-state wave functions 
used in the CSM. One also sees, in going from oxygen to calcium to lead at intermediate- 
to-high values of q, that the agreement between the two models becomes better and better. 
The slight shift in peak position seen in Fig. 3 for oxygen is seen to go away for the heavier 
nuclei and the importance of the tails of the response functions produced by the CSM 
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decreases with increasing A. With regard to the latter, in effect, the heavier the nucleus 
the more it appears to be dominated by volume effects, rather than surface effects, for the 
inclusive responses. At q = 300 or 500 MeV/c and low energy we see somewhat irregular 
behaviour in the shell model responses which is a consequence of the contribution of the 
potential resonances in that region. This behaviour is more pronounced in the case of 40 Ca 
than for 16 0, and, in the case of 208 Pb, it is responsible for the peak that clearly can be 
seen at q = 300 MeV/c and u = 40 MeV in both R L and R T . 

3. The rather good agreement between the RFG and CSM responses at high q is only obtained 
when the approximations discussed above are made. If the currents and kinematics are 
not "relativized" as we have in this work, but instead are taken to be their usual non- 
relativistic forms, then significant differences between the RFG and (conventional) CSM 
will be seen to occur. 

4. When comparing the five different electroweak responses, the worst agreement between 
the models is seen for R AV , although even there it is reasonably good for high q. It should 
be noted that this response is strongly suppressed in such mean-field calculations for the 
reasons discussed in Refs. || [|; indeed, at q — 500 MeV/c for example, it is only a few 
percent of the leading PV response R T AV (and the other, Ry A) is roughly 20% of this). As 
a consequence any differences seen here are unimportant for the PV asymmetry. On the 
other hand, especially at relatively low q this is known to be only part of the full story. 
In particular, it has been argued |13| that correlation effects that fall outside the scope 



of the strict RFG and CSM approaches can strongly influence this particular response at 
low-to-moderate momentum transfers, at the same time only mildly effecting the other 
four response functions. In the next subsection we return to this issue when discussing 
the role played by the electric strangeness content of the nucleon. 

5. In the figures we have shown the RFG results otained with the different choices for the 
Fermi momenta discussed above. The response functions are seen to differ somewhat, 
with a slight tendency to favor smaller values for pp. As we shall see immediately below, 
even this mild dependence on the Fermi momenta is considerably lessened for the PV 
asymmetry. 

Next we turn to the PV asymmetry. As seen in Fig. 9, such ratios are very insensitive to the 
choice of nuclear model made for the dynamics. Only in kinematic regions that correspond to the 
tails of the response functions do we see significant differences — and there the figure-of-merit 
for measuring the PV asymmetry is very small (see Ref. || for a discussion of the figure-of- 
merit). As a consequence, we show only the regions where the figure-of-merit is appreciable in 
the figure. 



17 



In the next subsection we shall discuss the sensitivity of PV observables to specific aspects 
of the nucleon's form factors. Following Ref. |H| we shall do this using the sum-rule asymmetry 
(SRA) defined as: 

^M^^"'^, (87) 



where the functions X' T and X' T are given in Ref. JT3[. These functions depend on the Fermi 
momentum and for simplicity we have used only the values (in MeV/c) pp = 225 ( 16 0), 235 
( 40 Ca) and 250/288 (p/n 208 Pb). The dependence on pp here is negligible and so choosing any 
reasonable value is an excellent approximation. The motivation for using the specific derived 
quantity 1Z is that it has been designed to be effective in removing most of the rapid variation of 
the single-nucleon form factors, to the extent that this is possible. While constructed with the 
RFG in mind, however, it should be noted that it merely leads to a reduced quantity in much the 
same way that dividing the nuclear cross section by the Mott cross section leads to a form factor 
- there is no loss of information, only a re-expression of the same information. Experimentally 
this quantity can be constructed from the measured helicity-difference and helicity-sum cross 
sections (PV and EM, respectively) upon dividing by the X's and integrating over a range in u. 
For the RFG the range of integration is naturally taken to include the entire response, since the 
sharp cutoff guarantees specific limits. In the case of the CSM, we integrate from u at threshold 
up to sufficiently large values of uo that the integrals in the numerator and denominator of 
Eq. ( ^7[ ) saturate. 

In Fig. 10 we show the SRA for 40 Ca and 208 Pb (the 16 results are very similar to the 
former and so are not shown). The sets of curves shown are described in the figure captions (see 
Fig. 4 for the basic key). Rather clearly the nuclear model dependence is very weak, especially 
for intermediate-to-high values of q, in that the various curves cannot be distinguished in the 
figure. This weak nuclear model dependence is essential if one wishes to use nuclei to explore 
properties of the nucleon itself. 

In Fig. 11 we show the CSM predictions for the SRA now for all three nuclei, observing 
the near universality of this observable for a large range of nuclei with only slight differences 
appearing when lead is compared with oxygen and calcium (the last two are essentially indis- 
tinguishable). 



3.3 Sensitivity to variations in the form factors of the nucleon 

We now turn to an assessment of the nuclear model dependences seen using the CSM and RFG 
when attempting to extract information about the isovector axial-vector and electric strangeness 
form factors of the nucleon. We do this using the sum-rule asymmetry defined above. 

We begin by presenting results in Fig. 12 for the SRA for 40 Ca and 208 Pb with no magnetic 
strangeness (as in the previous figures) and with fj, s in Eq. fl38|) set equal to —0.35. The latter 
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is chosen in accordance with the discussions of modeling the nucleon's strangeness (see || 
and references therein). For comparison, as also discussed in that reference, we note that the 
SAMPLE experiment being performed at Bates to measure backward-angle PV ep scattering 
is expected to have an uncertainty of approximately \Afi s \ ~ 0.2. In the figure we see that 
the nuclear model dependence is weaker than the dependence on G^} for intermediate-to-high 
values of q in the case of an N = Z nucleus such as 40 Ca ( 16 is almost indistinguishable from 
40 Ca). However, at low q or for nuclei with iV significantly larger than Z, such as 208 Pb, the 
dependence is much weaker. As discussed in Ref. 0, this arises from the fact that in the 
PV asymmetry this form factor is multiplied by the combination [ZG V M + NGm)/A, which is 
proportional to (Zp p + Np n )/A = 0.44 (—0.06) for calcium (lead). 

For backward-angle PV QE electron scattering the asymmetry depends importantly on the 

(s) 

magnetic strangeness form factor G M discussed above and the isovector, axial-vector form factor 
G^~^ parametrized in Eq. ([36]). The electric strangeness form factor G^ is unimportant for 6 e 

(T=l) 

large. In Fig. 13 we show the dependence on G A . The three sets of curves shown correspond 
to results with the isovector, axial-vector coupling at Q 2 = taken to be its canonical value, 
Qa = 1.262, together with results when qa = 1.262±10% (indicated with a (±) in the figure). In 
all cases the strangeness content of the nucleon has been neglected. Clearly, as in Ref. ]13[ where 
other nuclear models were inter-compared in the same way, for intermediate values of momentum 
transfer the nuclear model dependences are quite weak compared with the level of variation in 
1Z caused by 10% changes in the isovector axial-vector coupling. Hence the conclusion found in 
previous work still stands that measurements of this observable could constrain the NC axial- 
vector form factor of the nucleon with interesting precision. It is also clear from Figs. 12 and 
13 that measurements on two nuclei are needed to determine both fi s and g A — PV elastic 
scattering from the proton and PV QE scattering from nuclei such as those discussed in the 



present work or the deuteron discussed in Ref. |28] would provide the required information. 

Turning finally to a discussion of forward-angle PV QE electron scattering we focus on the 
dependences in the SRA on G'^, G$ an d G^ (G r ^ _1 ' ) is unimportant for 6 e small). In Fig. 14 
TZ is shown for three values of the Galster parameter X n in Eq. fl34"|), corresponding to ±10% 
variation in G% at its peak value. This sets a scale for the uncertainty expected in the SRA 
at forward angles, being typical of the hoped-for precision in determining G% from polarization 
measurements of parity-conserving e 2 H and e 3 He scattering. The Gj$ sensitivity has been 
shown in Fig. 12 and the dependence on G^ is shown in Fig. 15. In the latter three models 
are employed: (I) no electric strangeness; (II) G$ with p s = -3, = 5.6 in Eq. ©; and 
(III) Gg with p s = —3, A# = 0. Model (II) corresponds to "modest electric strangeness", 
while model (III) has "large electric strangeness" (see Ref. ||). Clearly when the interest is 
to see effects from electric strangeness the momentum transfer cannot be too low, since the 
single- nucleon form factor G% vanishes as |Q 2 | — > (see Eq. (|37|)). In the figure we have chosen 
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to present results for 40 Ca and 208 Pb at q = 300, 500 and 700 MeV/c; again the 16 and 40 Ca 
results are rather similar and so the former are not displayed. We see that at forward angles 
and high momentum transfer the sum-rule weighted asymmetry shows clear signatures of the 
electric strangeness, while at lower values of q the nuclear model dependences become important 
enough for the case of 208 Pb to obscure the G$ effects. Comparing Figs. 12, 14 and 15 it is 
clear that, given PV ep and QE information at backward angles to determine G^ and G^~ 1 ^ 
as well as possible, forward-angle PV QE scattering at high q should be capable of determining 
Gg\ For low q the situation is much less clear. For 40 Ca the results in Fig. 15 would appear 
to indicate at least that models II and III could be distinguished from model I for the electric 



strangeness. However, it is already clear from the results presented in Ref. |13|] that correlations 
beyond the mean field play a dominant role at low q. In recent work [2{J these effects have been 
explored in depth and, while the determining factor in the SRA at forward angles and low-g, 
they become sufficiently weak at intermediate-to-high values of q that the above conclusions 
concerning the dominant sensitivity to remain valid. 



4 Conclusions 

In conclusion, in this work we have (1) developed the single-nucleon vector and axial- vector 
current operators for use in parity-conserving and -violating electron scattering studies, (2) 
explored approximations to these currents and to specific aspects of relativistic kinematics within 
the context of the Fermi gas model, and (3) applied the same ideas to studies of electroweak 
QE responses and asymmetries in the continuum shell model. 

In developing the current operators we have seen that the traditional expansions in all 
dimensionless momenta provide poor approximations to the correct on-shell answers at high 
momentum transfers. By making more limited expansions in rj = p/M, but not in k = q/2M, 
A = uj/2M or t = k 2 — A 2 we have obtained approximate currents that should be appropriate in 
the region where the QE responses are large. Naturally, if r\ is of order unity or large, then this 
scheme will be suspect; however, in such circumstances the idea of using any on-shell current 
must also be questioned and we must resort to off-shell prescriptions for the currents. 

To evaluate the quality of the approximations made to the currents we have considered 
the relativistic Fermi gas model where exact (model) results exist. We have verified that the 
approximate currents are very successful for the closed shell (spin-saturated) nuclei studied in 
this work. Furthermore, we have verified that the non- relativistic Fermi gas model extended by 
incorporating the approximate currents and as well by "relativizing" the kinematics provides 
excellent agreement with the RFG results for the electroweak responses. This has led us to 
propose that the same approximate currents and assumption about the kinematics be invoked 
for any non-relativistic model of the nuclear dynamics involved in QE electron scattering. 
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In particular, we have made these assumptions in computing the QE responses with the 
continuum shell model. The results are very encouraging: the resulting electroweak response 
functions are rather similar to those found with the RFG model when the momentum transfer is 
high enough for the latter to be valid. Said the opposite way, if the two assumptions concerning 
the currents and kinematics were not to be made then the CSM and RFG results would be 
dramatically different at high q. 

Finally, we have employed the RFG (with various choices of Fermi momenta) and the CSM to 
study the PV asymmetry. We find very little nuclear model dependence in the results, especially 
those relating to the sum-rule asymmetry. This confirms the belief from other work that at high 
q the SRA can be used to study specific properties of the nucleon itself without incurring too 
much uncertainty from the nuclear modeling. Backward-angle PV QE scattering together with 
PV elastic ep scattering will serve to determine the magnetic strangeness and isovector, axial- 
vector form factors of the nucleon, while forward-angle PV QE scattering at high-g may add 
information about the electric strangeness form factor of the nucleon. 

Appendix A: On-Shell Single-Nucleon Electroweak Current Operators 

The on-shell single-nucleon four-vector EM current may be written in the form: 



J"(PA;P'A') =u{P'A' 



u(PA), (58) 



where the incident nucleon has four-momentum P^ = (E,p), the outgoing nucleon has four- 
momentum P' M = (E',p') and Q 11 = P' M — P M . The spin projections are labeled A and A' 



for incoming and outgoing nucleons, respectively. As in previous work fL0| , the dimensionless 
variables 

A = tojlM 
k = q/2M 

r = k 2 -X 2 (59) 
77 = p/M 

n 2 



e = E/M = ^/l + r] 

prove to be convenient. Let us also introduce the angle 9 between k, and r\. We then have the 
following relationships amongst the kinematic variables: 

kt) cos 9 = Xe — t 
T (e + Xf = k 2 (i + t + 5 2 ), (60) 
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where we define 5 = r/smO. For reasons that will become apparent later we also define the 
following quantities: 



A*2 = 



K 



vTT 



r(e + X) ^/TT^p 
2/?v / TTr 



1 ^ k r"l 



(61) 
(62) 



To evaluate the current in Eq. (pq) we insert the appropriate 7-matrices between the spinors 



u(P'A') 
u(PA) 



y/2 
1 

71 



1 + e 



7*7 



vTT 



Xa, 



(63) 



V Thv-o- J 



where Xa an d Xa' are the usual 2-component spin-| spinors, with rj' and e' the final outgoing 
nucleon analogs of the unprimed quantities. We wish to have expressions for the single-nucleon 
EM current operators J M (P; P') that occur inside these latter quantities, viz. 



J»(PA;P'A r )= ) &J»{P;P r )XA- 



(64) 



We use the bar over the current in the present appendix to distinguish an operator from its 
spin matrix elements; in the body of the paper for simplicity we suppress this notation. Writing 
these in the following way with an overall factor / removed (note that is not a fo ur- vector) , 



3» EE f V» 

fo = 



^V 1 + 4{I^2 2 6 2 

the EM current operator may then be expressed in terms of the quantities defined above, 



(65) 
(66) 



V° = £o + i£o(*x»7)-<T 

V 3 = (X/k)V° 



v 



6 



T] 



K ■ T] 



(k x tr) 



(67) 
(68) 
(69) 



77 



K ■ T] 



+£> ( K • <r) (k x 77) + [(« x 77) ■ cr] 

where the £'s (no spin dependence) and £"s (spin dependence) are the following: 

G E + 



^ = -r 

\ T 



-6 tG m 



2(1 + r 
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Ml^fi + ^2tGm 



« V 2(1 + r 

V 



-<r G 



(70) 



2k 2 v / TT7 



r/^l + 4/xi(/x 2 - Mi) - (niUzS) G 



M 



2k 3 



■H1H2G 



M 



2k(1 + r 



-/ii/i 2 [C^ — G 



Af I 



These expressions constitute exact expressions for the on-shell EM current operator. Equation 
( |S8f) reflects the fact that the current is conserved. 

The vector NC analogs of the results here are simply obtained by replacing Ge by Ge and 
Gm by Gm- 

The NC axial- vector current may be written as in Eq. (j58|) : 



4(PA;P'A') = u(P'A' 



j 5 u(PA), 



(71) 



where the "5" is used to indicate axial-vector quantities (from the extra 75 above, compared 
with the EM current in Eq. (|58|) ) and where Ga and Gp are the axial and induced pseudoscalar 
form factors of the nucleon, respectively. The analog of Eq. (E4l) becomes 



(72) 



with 



> 5 = foA", (73) 

where /o is given in Eq. (|66|) (as above, here A^ 1 is not a four- vector) . Clearly the pseudoscalar 
contributions are absent for transverse (1 and 2 components) projections of the current, since the 
coordinate system used has been chosen to have the momentum transfer along the 3-direction. 
One finds from carrying out the procedures discussed above for the EM (vector) current, inserting 
the spinors and 7-matrices, that the transverse projections are given by 







77- 









<i (kxi?) G a 



(74) 



with 



Ci 
CI 



2 "2"U, 

\/l +T/J.I, 



MlM2 



(75) 
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in parallel with Eqs. flTop . 

The exact results given above in Eqs. (|64] ) -(|7"0| ) may be used in treating EM interactions with 
an on-shell nucleon. However, for applications to nuclear physics as in the present work where 
nucleons are off-shell in general it is necessary to make some approximations. What is often done 
is to continue to use the above expressions as operators inserted between nuclear wave functions, 
where the latter are not on-shell plane waves but are the single-particle wave functions for 
nucleons in nuclei (i.e., in the presence of the nuclear mean field). The coordinate space operator 
corresponding to r\ is then — tjV. For typical situations the nucleon to which the virtual photon 
attaches is relatively low in energy — the dimensionless momentum 77 ~ rjp = Pf/M ~ 0.25 
yielding dimensionless energy e = 1 + \rf ~ 1 + ^rjp ~ 1.03. Accordingly the usual practice is 
to expand the expressions for the EM current given above (at least) in powers of 77, retaining 
only the leading-order terms. Equations ([70],|75]) have been cast in forms that make this easy. 
Let us work to linear order in 77. First, from Eqs. ( plf ) and fl66|) we have \l\ = 1 + 0(r] 2 ) 



and f = 1 + Oijf ), whereas it may be shown that /i 2 = 1 + iyT^7 ? 7 cos ^ + ^(v 2 ) D Y using 
Eqs. (|60H62|). We then have for the EM current operators carried to order 77: 



J° = 


z—Ge + - 






J 3 = 


(A/k) J° 


J 1 = 









Gm — ~^Ge 



(k x 77) • <r + 0{rf) (76) 

(77) 



M 



[k x er] H — — — — - (k ■ <t) (k x t]) 



Ge + 2 T Gm 



2(1 + 7-) 



K ■ T] 

, n-\—^) K 



employing Eqs. (p7j-|69D and noting from Eq. (|60| ) that, after some work, k = yr(l + r) + 
T7] cos 0+O(rj 2 ). Of course, when computing matrix elements of these operators and then forming 
bilinear combinations of the results to obtain the EM observables terms of order rj 2 must be 
neglected if the operators themselves have been expanded only to order 77, since other terms will 
enter from considering the neglected 0(r] 2 ) contributions in Eqs. (|7"6"H78"D . 

Note that we have only expanded in powers of 77, having argued that this is typically a 
small quantity. Indeed, if it is not small (for instance, when considering very high momentum 
components in the nuclear wave function, viz. 77 ~ 1 «-> p ~ M) then the full expressions must 
be employed; however, in such circumstances the struck nucleon is very far off-shell and the 
entire procedure becomes doubtful. What we have not done is to expand in powers of the other 
dimensionless momenta, k, A or r. To see why this can be a problem at moderate momentum 
transfers note that the typical observables that occur (for example, W 2 = {G 2 E + tG 2 m )/{1 + r)) 
involve combinations of the form factors where one often has G\ together with {yfrG^if ' . Using 
the dipole approximation to estimate the size of these quantities we see that they become equal 
for protons at r S 0.13 \Q 2 \ = (670 MeV/c) 2 , for neutrons at r £ 0.27 «-> \Q 2 \ = (980 
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MeV/c) 2 , or alternatively, for isovector contributions at r = 0.045 <-+ \Q 2 \ = (400 MeV/c) 2 
and for isoscalar contributions at r = 1.29 <-> |Q 2 | = (2.1 GeV/c) 2 . Clearly for many practical 
applications it is not valid only to retain terms of leading order in ypr or k. Moreover, as we see 
from the results presented in Sect. 3.1, the distinction between k and yfr becomes important for 
QE scattering at high momentum transfer. In fact, as we see from the expressions given above, 
it is unnecessary to make expansions in k, A or r at all. If one does, however, choose to do so, 
then at intermediate momentum transfers the combination G' M = yfrGu should be regarded 
as being of leading order, and not of order M~ 2 as is often assumed. With these caveats under 
some circumstances Eq. (|78|) may be approximated by 



J L ^-^-{iG M [kxo-]-G e 

K 



k ■ rj 

7 '1 I ~^2~ J K 



(79) 



At this point it should be mentioned that neglecting the second term in Eq. (|76|) (the spin- 
orbit term), which is order nr] compared with the first term which is of order unity, may not 
always be justified. However, in spin-saturated systems such as the closed-shell nuclei considered 
in the present work this term contributes only quadratically and may safely be neglected when 
computing EM responses. For instance, the contribution of the spin-orbit term to R L and R\y 
in a spin-saturated system can be estimated using the RFG model to be ^rf F ^ (l + §f) and 

YoVft+^ (l + Ggi" )' res P ec ^i ve ly? compared with the leading-order contributions ^G 2 E and 
^tGeGe- For typical values of rj F (~ |) the correction to R L is negligible (< 4%), while the 
correction to R L AV is sizable (~ 40%) and has been taken into account in the results presented in 
Sect. 3. There we see that the latter response is quite small compared with "normal" response 
functions and thus this correction is relatively unimportant. 

As before, the vector NC analogs of these EM results are obtained simply by making the 
replacements Ge,m — ► Ge,m- Likewise, the axial-vector results above may be expanded in 
powers of rj following the procedures used for the EM case, yielding the transverse axial-vector 
current operator to order rj: 

Ji; = VT+^Ga^ 1 - 

1 



2(l+r) 











(W) 


rj - 




-*(« X j 







0) 



vrrfGA^, (si) 

where in Eq. the second term, involving corrections of order kt] compared with the first 
term, has been neglected. 

Finally, in the discussions in the main body of the paper we also refer to the strict non- 
relativistic limit. In the present work we take this to mean that the currents in Eqs. ( [76 , 79 , 81]) 
are futher approximated by taking k and r to be small compared with unity to yield 

J° — > G E (82) 
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J 
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_L 



-iG M [kx(t] + G e 
Ga^. 



k, ■ r\ 

K 2 



(83) 
(84) 



Appendix B: Reduced Matrix Elements in the Shell Model 

After a laborious but straightforward calculation involving angular momentum algebra, the 
different reduced matrix elements involved in the PV responses can be written as follows: for 
the EM matrix elements one has 

jh J 



(ph-\j\\M c J m \\0) = 



[ ^^^ G ^ + £h + J ^Hlj2 -1/2 



r-oo 

x / drr 2 jj(qr)R p (r)R h (r) 
Jo 



(85) 
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~ 1)Jp+1/2 ^G M i{i v + 4 + j)] P j h j f /; J * i 

4ttJ(J+ 1) V l / 2 ~ 1 / 2 



x 



ilp-lh) / drr 2 jj(qr)R p (r)R h (r) 
Jo 



(86) 



(^-viraio>c 



- (-i)^ v 2 ^ n , s f / j P j„ j 

yW(J+l) K V X / 2 _1 / 2 

X 2M i, drrjj ^q \ _ 7/l ^> + ^ + 1) + ^ + 1)] R P (r)-R h (r) 

(87) 



+ [(7 P - 7ft)(7p + 7fc + 1) " ^ + 1)] Rh(r)-^Rp(r) 



(p/i" 1 , J|E J ||0) mag = / ; VfG M £(£ p + £ h + J+ l) JpJh J 
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1/2 -1/2 



^/4ttJ(J+1) 

/•OO 

/ drr 2 jj(qr)- 
Jo i 



R P (r)R h (r) 



(ph- 1 , J||Tf/||0) conv = ' ( l)3P V2 ^G E Z(£ P + £ h + J + l)j p j h J 

y/AnJ(J+l) k 

l P /2 _ J * /2 q j [(7p + 7h)(7p + 7, + 1) - J( J + 1)] ^ /f drr^i^COi^Cr), 



(89) 
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together with their vector NC analogs, obtained by replacing Ge with Ge and Gm with Gm- 

For the axial- vector NC operators, using the relations given by Eqs. ([)4],|)5|), one simply 
obtains 



(ph-\J\\(T£ J )ji 

j P jh J 
1/2 -1/2 



'_l)i P -i/2 
AnJ(J + l] 



VI + tG a £(I p + 4 + J + l)j P JhJ 



oo "j 

dr r 2 jj(qr)- 
u q 



(7 P + 7.) ^; + -) + 



1\ J(J+1) 



Rp{r)R h {r) 



(90) 



/•oo 

x (7p~7fc) / drr 2 jj(qr)R p (r)R h (r). (91) 

w 

In all of the above equations we have used the sign function £ defined as 

£{1, + li + J) = • (92) 

The functions, R p (Rh) represent the radial single-particle (single-hole) wave functions, and we 
have introduced the term ji = — ji)(2ji + 1). 
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Appendix C: Fermi Gas Model 

In this appendix we summarize the main features of the relativistic Fermi gas model (RFG), 



drawing on the detailed discussions in Refs. ||, |10| [11], [12|, D|. The purpose in presenting this 



material again is twofold: first, we make comparisons with the RFG and its extensions |13| when 
discussing PV electron scattering in the results section of the present paper; secondly, there are 
approximations that we make when employing the CSM and the RFG approach provides the 
basic motivation for these, as discussed below. 

We begin by re-stating the procedures followed in deriving the RFG responses in the above- 
cited work. The exact, fully-relativistic, on-shell, single- nucleon current matrix elements ob- 
tained by inserting spin spinors before and after the current operators given in Eqs. (|65H 70l) 
and (|73|-|75l) may be computed in a straightforward way. The results are given in Ref. |Tl| (in 
fact, multiplied by factors that occur in the invariant phase-space factors). These are the matrix 
elements for elastic electroweak scattering from on-shell, moving nucleons. Taking, as usual, the 
RFG step-function momentum distribution characterized by dimensionless Fermi momentum 
rjp = pf/M and integrating over all nucleons in the filled Fermi sea, one obtains the familiar 
response functions 



R 
R 



L,T 



L,T 
AV 



R 



VA 



R (k, X)U l ' t (k,X) 
R (k, X)U l ' t (k,X) 

Ro(k, x)u t '(k, A), 



(93) 
(94) 
(95) 



where the dimensionless momentum and energy transfers are defined in Appendix A. One may 
consider both the Pauli-blocked and non-Pauli-blocked regimes (k < r\p and k > tjf, respec- 
tively). Focusing here on the latter and following Ref. [HJ we define a scaling variable 




(96) 



where £p = J 1 + rjp — 1 = r/p/2, since rjp is typically small, as discussed in Appendix A. The 
RFG QE peak occurs at A = A = \ Wl + 4k 2 - 1]. 

In the non-Pauli-blocked regime the overall responses are characterized by 



Ro(k, A) 



^ 2 )^(1-^ 2 



(97) 



where Af is Z (N) when protons (neutrons) are involved (the expressions for the Pauli-blocked 
regime are discussed in Ref. [fllj - here our main emphasis is on the non-Pauli-blocked region 
although in Sect. 3 some results are also presented for k < rjp)- The dependences on the 



28 



single-nucleon content in the problem are then isolated in the factors 

.2 



U l (k,X) 

U t {k, A) 

U l (k, A) 
r ^,A) 



- {[(1 + r)W 2 {r) - W 1 (t)} + W 2 {r)A} 

T 

W 2 (t)A} 



WAt) 
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U t '(k,X) 



-{\(l + r)W 2 (r 

2Wi(r) + W 2 (r)A 
2W 1 (r) + W 2 (r)A 
t(1 + t)W 3 (t){i+A}. 



(98) 

(99) 

(100) 
(101) 
(102) 



Here the nucleon form factors are those discussed in the main body of the text (for protons or 
neutrons, as is appropriate). The factors A and A were introduced in Refs. [|, [TCfl : 



A 



A 




(103) 



(104) 



Clearly these two terms provide relatively small contributions to the totals in Eqs. (|98Hl02j ), 
since they both involve the factor £p which is very small. These are contributions of order if F 
and typically enter at the few percent level in the kinematic regime of interest here. Indeed, in 
Appendix A we argued that such 0(rj 2 ) effects could be neglected — we have done so in the 
treatment of the CSM in the present work. In Sect. 3.1 we present results for the RFG (and these 
agree with the past results cited above) both with the exact relativistic current matrix elements 
and with the approximations discussed in Appendix A. Moreover, we also display in Sect. 3.1 
some results for what we defined in Appendix A to be the strict non-relativistic approximation 



(see Eqs. (^2-84)) to illustrate how poor those commonly-made approximations can be (see 
Fig. 2). 

A second point to be drawn from this summary of the RFG formalism, in addition to bringing 
out the roles played by the single-nucleon current matrix elements, is seen by returning to Rq 
defined in Eq. ( |97|) . The basic RFG behavior in the non-Pauli-blocked region when the single- 
nucleon factors (the IPs and IPs) are removed is parabolic in the scaling variable ip. For the 
RFG this quantity ranges between —1 and +1, with the QE peak occurring when ip = 0. An 



excellent approximation to the scaling variable may be shown to be (see Ref. fl~2|l ) 



Vf 



A(l + A) 



K 



(105) 



making it easy to see that the QE peak (ip = 0) corresponds to A(l + \)=h?^\ = t^uj 
\Q 2 \/2M, the relativistic kinematical relation. The non-relativistic version of Eq. ( |105| ) is 



Vf 



K 



(106) 
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and yields the non-relativistic relation A = k 2 <-> oo = q 2 /2M for the peak position. In the 
past Fermi gas work cited above it was shown that a very good approximation is to employ 
the non-relativistic many-body responses (of course, using good single-nucleon currents), but 
with the replacement A — > A(l + A) or, for dimensionful variables, 00 — > uo(l + w/2M). The 
replacement is made only in the many-body aspects of the problem, not in the single-nucleon 
form factors which are evaluated at the correct value of r, namely that given by the electron 
kinematic variables. In Sect. 3.1 we present results to illustrate the quality of this replacement 
in the nuclear kinematics. With this as motivation, we follow the same "relativizing" procedure 
for the kinematics of the CSM as well by using the relativistic relation in the energy-conserving 
5-function by making the same replacement. 
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V [MeV) V ls [MeV] R [fm] a [fm] R ls {fm] a la [fm] 



16 p 


52.5 


7.00 


3.20 


0.53 


3.20 


0.53 


11 


52.5 


6.54 


3.20 


0.53 


3.20 


0.53 


40 Ca p 


57.5 


11.11 


4.10 


0.53 


4.10 


0.53 


n 


55.0 


8.50 


4.10 


0.53 


4.10 


0.53 


208p b p 


60.4 


6.75 


7.46 


0.79 


7.20 


0.59 


11 


44.3 


6.08 


7.46 


0.66 


6.96 


0.64 



Table 1: Parameters of the Woods-Saxon potential. They have been taken from Refs. [JIB] for 
16 0, 121 for 40 Ca, and IE3I for 208 Pb. 



Nucleus 


g[MeV/c] 


J max 


J con 


16Q 


300 


15 


8 




500 


15 


12 




700 


20 


20 


40 Ca 


300 


15 


11 




500 


15 


15 




700 


23 


23 


208 pk 


300 


20 


18 




500 


23 


23 




700 


30 


30 



Table 2: Multipoles used in the calculation. 
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FIGURE CAPTIONS 



1. One boson exchange diagrams considered in this work. 

2. Electromagnetic response functions R L and R T for 40 Ca shown at q — 500, 700 and 1000 
MeV/c: (solid curves) full RFG model; using the approximate currents given in Sect. 2.2, 
(long-dashed curves) with the replacement A — > A(l + A) made to the non-relativistic FG 
model as discussed in Appendix C or (dot-dashed curves) without this replacement; as in 
the previous case, but using the non-relativistic currents obtained by setting the factor 
k/^/t — > 1 in Eqs. (|20|) and (|22|) with (short-dashed curves) or without (dotted curves) 
the A — > A(l + A) replacement. In all cases pf is fixed at 235 MeV/c. 

3. Electromagnetic response functions R L and R T for 16 at q = 300, 500 and 700 MeV/c: 
(solid curves) CSM; (other curves) full RFG model, (dashed curves, pf = 215 MeV/c) and 
(dotted curves, pf = 225 MeV/c). 

4. Electromagnetic response function R L for 40 Ca and 208 Pb at q = 300, 500 and 700 MeV/c. 
The labeling in this and most of the following figures is: (solid curves) CSM; (other curves) 
full RFG model. The latter are the following: for 40 Ca two are shown, (dashed curves, 
Pf = 235 MeV/c) and (dotted curves, 245 MeV/c); for 208 Pb three are shown, (long-dashed 
curves, pp v = 250 MeV/c, pp n = 288 MeV/c), (short-dashed curves, pp p = PFn = 250 
MeV/c) and (dotted curves, pf p = 230 MeV/c, pf u = 265 MeV/c), as discussed in the 
text. In fact, for the longitudinal response the short- and long-dashed curves essentially 
coincide. 

5. As in Fig. 4, except now for the EM response R T . 

6. As in Fig. 4, except now for the PV response R\y- 

7. As in Fig. 4, except now for the PV response R v \. 

8. As in Fig. 4, except now for the PV response R\y. 

9. Parity-violating asymmetry for 16 0, 40 Ca and 208 Pb at q = 300, 500 and 700 MeV/c and 
at 6 e = 10° and 150° (labeling of curves as in Fig. 4). The difficulty in distinguishing the 
curves shows the nuclear model insensitivity of the results. 

10. The sum-rule asymmetry (SRA) defined in Eq. fl57j) shown for 40 Ca and 208 Pb at q = 300, 
500 and 700 MeV/c as a function of scattering angle 9 e (labeling of curves as in Fig. 4). 
The difficulty in distinguishing the curves, especially in the case of 40 Ca, reflects the high 
degree of nuclear model insensitivity for this observable. 
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11. The sum-rule asymmetry (SRA) defined in Eq. (p7| ) shown for 16 (solid curves), 40 Ca 
(dashed curves) and 208 Pb (dotted curves) at q — 300, 500 and 700 MeV/c as a function 
of scattering angle 9 e . The CSM is used in all three cases. The difficulty in distinguishing 
the curves (the results for 16 and 40 Ca are essentially equal) indicates the universality of 
this observable for a wide range of nuclei. 

12. As in Fig. 10, but now showing families of curves with magnetic strangeness (with fi s = 
—0.35) and without magnetic strangeness (indicated with a "0" in the case of 40 Ca where 
the families of curves can be distinguished at high q) . In this case the curves are seen to 
coalesce into sets determined by the value of fj, s and not by the nuclear model variations, 
as discussed in more depth in the text. The families of curves are labeled as in Fig. 4. 

13. As in Fig. 10, but now focusing on scattering at backward angles and showing families of 
curves with different values of the isovector, axial- vector form factor of the nucleon. The 
curves marked with a (±) correspond to setting = 1.262±10%, while the remaining 
curves have this quantity set to 1.262. The families of curves are labeled as in Fig. 4. At 
intermediate-to-high values of q the curves are seen to group into sets determined by the 
axial-vector coupling and not so much by the nuclear model variations. 

14. As in Fig. 10, but now focusing on scattering at forward angles and showing families 
of curves with different values of the neutron electric form factor, the standard Galster 
parameterization with A n = 5.6 (see Ref. ||) and values of this parameter that yield ± 
10% variation in G\ at kinematics where the form factor peaks, namely A n = 4.2 and 7.3 
(marked (±) in the figure when the groups are distinguishable). The families of curves are 
labeled as in Fig. 4. The nuclear model dependences and the differences due to variations 
in G% at this level would clearly be quite difficult to separate. 

15. As in Fig. 14, but now showing families of curves with different parameterizations for 
G [ e ] : the curves marked (I) have no contributions, marked (II) have p s = —3 and 
A^ = 5.6 ("modest electric strangeness"), while those marked (III) have p s = —3 and 
\g = ("large electric strangeness"). In all cases the magnetic strangeness form factor 
has been taken to be the standard parameterization discussed in the text. Clearly for high 
q and forward angles the nuclear model dependences are considerably weaker than the 
dependence on the parameterization of electric strangeness in the nucleon. 
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